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Abstract 

We study the effective resistance of small-world resistor networks. Utilizing recent 
analytic results for the propagator of the Edwards- Wilkinson process on small-world 
networks, we obtain the asymptotic behavior of the disorder-averaged two-point 
resistance in the large system-size limit. We find that the small- world structure 
suppresses large network resistances: both the average resistance and its standard 
deviation approaches a finite value in the large system-size limit for any non-zero 
density of random links. We also consider a scenario where the link conductance 
decays as a power of the length of the random links, In this case we find that 
the average effective system resistance diverges for any non-zero value of a. 
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Resistor networks have been widely studied since the 70 's as models for con- 
ductivity problems and classical transport in disordered media [1,2,3,4]. Re- 
lated studies on fuse networks have been investigated on random percolating 
lattices with various applications to breakdown processes in condensed matter 
and materials systems, ranging from brittle fracture to dielectric breakdown 
[5,6,7,8,9]. 
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Recent research on complex networks [10,11,12] has turned to focus on dy- 
namics on networks with apphcations to synchronization in natural and arti- 
ficial systems [13,14,15,16,17,18], and transport phenomena [19,20,21,22]. In- 
teresting recent studies have examined the tradeoffs between redundancy and 
pleiotropy [23], and centralized versus decentralized design [24], in complex 
networks. Finding the resistance between any two points on a complex network 
is tractable and builds upon early mesh-rcsistance techniques [25]. Estimating 
the strength of collaborative ties between nodes in collaboration networks [26] 
and quantifying the centrality of a node in weighted networks can also be mod- 
eled by resistor networks [27] . While resistor networks have been employed to 
study and explore community structures in social networks [27,28,29], they 
have not been investigated as prototypical models for transport phenomena 
in complex networks until very recently [30,31]. The work by Lopez ct al. [30] 
revealed that in scale-free (SF) networks [10,32] anomalous transport proper- 
ties can emerge, displayed by the power-law tail of distribution of the network 
conductance. 

Here we investigate the effective system resistance of small-world (SW) net- 
works [33,34,35]. Our results, in part, are based on recent calculations [36,37,38,39] 
of the disordered averaged propagator of the Edwards- Wilkinson (EW) [40] 
process extended to a SW network. The EW process on a network can be 
thought of in terms of a synchronization paradigm in a noisy environment. As a 
linear approximation, it also serves as the simplest model for generic causally- 
constrained queuing networks [41], such as manufacturing supply chains, e- 
commerce based services facilitated by interconnected servers [42] , and certain 
distributed parallel schemes on computer networks [18,43]. In the context of 
the latter, it was shown [18] that when extending the original short-range con- 
nections to a SW-like network (essentially, by adding a small density of random 
links on top of a regular graph), the spread between completion times of tasks 
performed on different nodes of a computer network remains bounded, rather 
than diverging over time. Further, an infinitesimal extra "cost" is sufficient to 
achieve this reduction. An important measure of efficiency is the spread (or 
"width") of task-completion landscapes in such processing networks (larger 
spread corresponds to longer delays and poorer efficiency). It is evident that 
this measure — the width of the EW landscape on a network — is identical to 
the average resistance (characterizing transport efficiency) of the same net- 
work. While this connection between the network propagator and the network 
resistance [2,44,45], just like the one between random walks and network re- 
sistance [46,47,48], is well known, it has not been exploited to study transport 
efficiency of SW networks. Further, the connection between the average spread 
of an EW steady-state landscape and the resistance of the same network gives 
some insight in treating synchronization and transport efficiency on the same 
footing. Namely, understanding the effects of the SW links in suppressing the 
diverging long-wavelength modes of the network propagator, originally present 
in regular lattices. 
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Our main result is that in SW networks, the average system resistance be- 
comes finite for an arbitrarily small density of random links, governed by the 
same behavior of the network propagator which is responsible for suppressing 
"rough" synchronization landscapes [18,36]. 

The Edwards- Wilkinson process on a network. — The EW process in a synchro- 
nization context on a network, is given by the Langevin equation 

TV 

9thi = -J^'^ijihi - hj) r]i{t) , (1) 



where hi{t) is the general stochastic field variable on a node (such as fiuctua- 
tions in the task-completion landscape in certain distributed parallel schemes 
on computer networks [18,36]) and rii{t) is a delta-correlated noise with zero 
mean and variance {r]i{t)r]j{t'))=26ij6(t — t'). Here, Aij=Aji>0 is the effec- 
tive coupling between the nodes (y4jj=0). Defining the network Laplacian, 
Fij — 5ij J2i -^ii ~ -^ij, we can rewrite Eq. (1) 

N 

dthi^-J2^ijhj + r]i{t) . (2) 
For the steady-state equal-time two-point correlation function one finds 

N-l 

Gij = {{hi - h){hj - h)) = fr.i = J2 T-'ipkiipkj , (3) 

fe=l '^k 



where h — (l/N) J^iLi K and (...) denotes an ensemble average over the noise 
in Eq. (2). Here, T^^ denotes the inverse of F in the space orthogonal to 
the zero mode. Also, {ipki}j^i and Afc, k = 0,1, . . . , N — 1, denote the A;th 
normalized eigenvectors and the corresponding eigenvalues, respectively. The 
k — index is reserved for the zero mode of the Laplacian on any network: 
all components are identical of this eigenvector and Aq = 0. The last form 
in Eq. (3) (the spectral decomposition of F"^) is useful for exact numerical 
diagonalization purposes. As one can see from Eq. (3), G is the inverse of the 
coupling matrix F in the space orthogonal to the zero mode of the Laplacian. 
In particular, the average spread or width in the synchronization landscape 
becomes 

/ . N \ . N 1 1 

(.^).(-g(..-.)^).-gC..^-g3^. (4) 



For large networked systems, the above observable is typically self-averaging 
(w^) ~ [(w^)], where [. . .] denotes the average over the network disorder. Thus, 
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if one is able to calculate the disorder- averaged propagator [Gij] , it provides 
the scaling behavior of the average spread of the synchronization landscape in 
the limit of iV— >oo. 

The two-point resistance of a network. — The stationary currents and voltages 
in any network of resistors are governed by Kirchhoff 's and Ohm's laws 

n 

where is the conductance of the link between node m and n, and 1^, is 
the net current flowing into the network at node m. Note that is zero, 
unless node m is connected to an external terminal. Connecting the network 
to a "battery" with fixed voltage drop V through nodes i and j as the input 
and output terminals, yields 

X] TmnK = I{5mi " ^mj) : (6) 
n 

where Tmn is the same network Laplacian as introduced earlier in the context 
of the EW process [Eq. (2)], associating the link conductance with the coupling 
matrix there; 5ij is the Kronecker delta. Here, / is the magnitude of the current 
entering and leaving the system at node i and node j, respectively. Solving for 

the voltages is well defined in the subspacc orthogonal to the zero-mode of the 
network Laplacian; the right-hand side vector of Eq. (6) is in that subspace. 
Hence, introducing the voltages measured from the mean = Vjn — V, where 
V={l/N) J2m=i ^rni and employing F"^ one has 

n n 

where G is the same network propagator discussed in the previous section 
in the context of the EW process on networks [Eq. (3)]. Applying the above 
equation to the voltage difference across node i and node j to which the battery 
is attached, one finds 

V ^Y,-Y^ = %- Vj = I {Gu + Gjj - 2Gij) . (8) 

For the equivalent two-point resistance between node i and j one finally ob- 
tains 
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where the last form in Eq. (9) is, again, useful for exact numerical diago- 
nalization purposes. Looking at Eq. (9), one can realize that the two-point 
resistance of a network between node i and j is the same as the steady-state 
height- difference correlation function of the EW process on the network 

{{hi - hjf) = Gii + Gjj - 2Gij = Rij . (10) 

The height-difference correlation function is a standard observable in surface- 
growth phenomena, extensively studied in the past two decades [49], so many 
of the answers for regular resistor networks can be obtained directly by looking 
at the equivalent EW model on a rf-dimensional substrate. For example, in an 
infinite one-dimensional system, the resistance between two nodes, separated 
by a distance |i — i|, diverges with the separation as Rij = R{\i — j\) ~ \i — j\ 
[44,45,49]. Another trivial, but insightful, relationship between the EW process 
and the resistor network can be obtained by summing up Eq. (10) over all i ^ j 
pairs, yielding 

^-]v(Ari^g^--2(^'>' (11) 

i.e., the average system resistance of a given network is twice the steady-state 
width of the EW process on the same network. 

Effective Resistance of simple SW Networks. — First, we consider "simple" SW 
resistor networks, where the conductance of each link is identical, with unit 
value, for simplicity. When studying network-transport phenomena for systems 
where physical links are subject to strong cost and geometric constraints, this 
can be unrealistic and cost- prohibitive. For others, e.g., modeling information 
flow in social networks [27,28,29,50], this can be an acceptable starting point, 
since "long-range" connections do not necessarily degrade the information car- 
rying capacity and the efficiency (e.g., influence) of that link. We start with 
a one-dimensional ring with nodes (i.e., impose periodic boundary condi- 
tions), and add a "random" link to each pair of nodes, independently for each 
pair, with probability p/N. In addition to the two nearest- neighbor connec- 
tions, now each node, on average has p random links, so p is the density of 
random links. The resulting network is essentially an Erdos-Renyi (ER) net- 
work [51] on top of a one-dimensional graph. This SW construction slightly 
differs from the original Watts-Strogatz one [33] where random links are in- 
troduced through "rewiring". The resulting network, however, has the same 
universal properties in the small-p, large- A" limit [52,53], and is also more 
amenable to analytic approximations. 

The coupling matrix for the differences of the relevant variables [Eq. (1) and 
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(5)] then becomes 

Aij = Sij-i + + Jij , (12) 



where the matrix elements Jij are quenched random variables; Jij=l with 
probability p/N and Jij=0 with probability (1 — p/N). The corresponding 
Laplacian then can be written as 

^ij = 2(5ij - Sij_i - Sij+i + Sij ^ Jii - Jij . (13) 

I 



Equations (12) and (13), with Jij defined above, correspond to identical (unit) 
conductance for each existing connection in the resistor network. Our numer- 
ical scheme relied on the exact numerical diagonalization of the SW network 
Laplacian F in Eq. (13) [54]. Our analytic results, asymptotically exact in the 
large system-size limit, are straightforward applications of those of the EW 
propagator on SW networks [36,37,38,39]. 

Averaging over the network- disorder restores translational invariance, hence 
the disorder-averaged two-point function [Gij] = [G{\i — will only depend 
on the underlying Euclidean distance between the nodes. These correlation 
functions have been calculated using disorder- averaged self-consistent pertur- 
bation theory [36,37,38,39]. For the disorder-averaged two-point function for 
small p values, in the infinite system-size limit one finds [36,37,38,39] 

m)] ^ ^e"^' > (14) 



where E ~ is an effective mass generated by the random links for simple 
SW networks [36,53]. Then, using Eq. (9), for the average resistance on a SW 
network between two nodes separated by a distance I, we obtain 

[R{1)] = 2 ([G(0)] - [cm ^ ^ (l - e-^') , (15) 



approaching a finite value in the limit of infinite separation, lim;_»oo[-R(0] — 
S^^/^ ~ p^^. In contrast, on a regular one- dimensional ring, the resistance 
between two nodes separated by a distance / diverges in a power-law fashion, 
R[l) ~ /, as can be seen by taking the S— >0 limit in Eq. (15) or by direct 
calculations on regular lattices [44,45] . Further, the average resistance is finite 
for an arbitrarily small but non-zero p in the limit of N^oo, 

Rc^[R]^ 2[{w')] = 2[G(0)] c^^^p-\ (16) 
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Fig. 1. Disorder-averaged (a) two-point function and (b) two-point resistance as a 
function of the separation / in simple SW networks for p = 0.10 and three sys- 
tem sizes. The soHd Hne in (a) and (b) corresponds to the exponential decay and 
saturation given by Eqs. (14) and (15), respectively 
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Fig. 2. (a) Average resistance vs the density of random links. The straight solid line 
indicate the asymptotic infinite system-size behavior [Eq. (16)]. (b) Average resis- 
tance vs the system size in simple SW networks. The straight solid line corresponds 
to the behavior of the one-dimensional regular network (ring), [R] ~ 

in strong contrast with average resistance for a regular network diverging as 
R ~ N/6. Equations (15) and (16) are the central results of our paper. They 
capture the average resistance between nodes separated by a distance / and the 
average system resistance for SW networks with a small but non-zero density 
of random links, respectively. Results from exact numerical diagonalizations, 
shown in Fig. 1 and Fig. 2, up to systematic finite-size effects, agree very well 
with the above predictions. 



In addition to the above asymptotic results, valid in the infinite system-size 
limit, we also constructed the scaling form, capturing the finite-size effects, 
e.g., for the average resistance [38,55]. From the above it is clear that in 
addition to the linear system size N, there is one other lengthscale in the 
problem for non-zero p values, ^ = l/y/T, ~ p~^. This lengthscale is, in fact. 
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Fig. 3. Scaled average resistance in simple SW networks according to Eq. (17). The 
straight solid line corresponds to the asymptotic behavior of the scaling function 
for small arguments [Eq. (18)]. 

the average distance between nodes which have random hnks emanating from 
them. For p = (the limit of a regular one- dimensional network) [R] ~ A^, 
while for p 7^ 0, in the infinite system-size limit, it approaches a constant, 
[R] ~ 1/VTi = ^ ~ [Fig. 2(a) and (b)]. Thus, the finite-size behavior of 
the average resistance can be expressed as 



[R] = Nf{^/N) , 

where f{x) is a scaling function such that 

. X if X <^ 1 
const, if X 1 . 



(17) 



;i8) 



The scaled numerical data, [R]/N vs ^/N [Fig. 3], shows good collapse, as 
suggested by Eq. (17). 

We also studied the probabihty distribution of the effective resistance of the 
network [Fig. 4]. The overall distribution is shown in Fig. 4(a). Further, we 
constructed the distribution of the effective network resistance between two 
nodes separated by a distance I, P{R\l) [4], indicating that they converge to 
a limit distribution for l—>-oo [Fig. 4(b) and (c)]. These results imply that for 
SW networks, both small and large effective resistance values are strongly (at 
least exponentially) suppressed about the average. This is in strong contrast 
with the behavior of SF resistor networks [30], where large resistance values 
are strongly suppressed, but the probability of small values decays only in 
a power-law fashion; hence a power-law tail in the conductance distribution 
occurs for large g=l/R values. This finding for SF networks implies [30], that 
there exist a few nodes ( "hubs" ) in the system that, if selected as the input and 
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Fig. 4. (a) Network resistance probability distributions (histograms) for p = 0.10 and 
for three system sizes, (b) Network resistance probabihty distributions (histograms) 
for nodes separated by a distance for various I values for N = 4000. (c) Behavior 
of the average [-R(/)] and the standard deviation [Ai?(/)] of the distributions shown 
in (b). 



output nodes, can support anomalously large transport through the network. 
This phenomenon is absent (as one can expect) in SW networks, just like in 
completely random (ER) networks [30] related to the exponential tail of the 
degree distributions of these networks. In Fig. 5 we compare the conductance 
distributions for the SW and the Barabasi-Albert (BA) [32,56] SF network, 
with the same average degree and uniform link conductance. At this point we 
note that while anomalously large conductances are absent in SW networks, 
they are more efficient, on average, in supporting transport between two arbi- 
trary pairs of nodes, i.e., [-R]5vk < [-R]sf, and \g\sw > [g]sF- In comparison, for 
the two networks shown in Fig. 5 with the same average degree and uniform 
link conductance, the average network resistance and conductance values are 
[R]sw ^ 0.472, [R]sF ^ 0.572, and [g]sw ^ 2.28, [g]sF ^ 1-93, respectively. In 
real-life complex networks with SF structure, however, the link conductances 
are typically weighted [21], ultimately leading to better performance for SF 
networks [30]. 

Effective Resistance of SW Networks with Distance-Dependent Conductances. — 
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Fig. 5. Comparison of the network conductance distribution with that of the BA 
scale- free network model with the same average degree {k)=6 (BA network with 
m = 3 [56], and SW network with p = 4.0.) 

Now we consider the more general case where the conductances of the random 
(possibly long-range) links decay with the underlying spatial distance between 
the nodes they connect in a power-law fashion: Jij = — j|" with prob- 
ability p/N and Jjj=0 with probability (1 — p/N) in Eq. (12). Keeping the 
density of the random links, p, fixed (at a non-zero value) and taking the 
large system-size limit, corresponds to the limit of fixed density of weak links 
at large scales. Then one can argue that, to leading order, mean-field scaling 
holds [37,38,39,55,57]. Focusing on the 0<a<l regime, we find that the aver- 
age link strength, decaying as [Jij] = {p/N)\i — j\~°', gives rise to a system-size 
dependent effective mass T^n — (-^/2)~"p/(l — a) and consequently 



Thus, the average system resistance diverges with the system size for an ar- 
bitrarily small but nonzero value of a. Figure 6 supports this picture, but 
also indicates that corrections beyond the mean-field approximation are im- 
portant and noticeable for the range of system-sizes that were accessible via 
numerical methods. In fact, an analysis of the naive perturbative approach 
[37] reveals that although higher-order corrections are becoming progressively 
smaller as increases, their prefactor is singular for certain values of a [55]. 
Given these subtleties, the deviations [Figure 6] from the predicted asymptotic 
scaling Eq. (19) are reasonable. 

For the interested reader, familiar with the diagrammatics of [36], we present 
a brief analysis of the higher order corrections. The mean-field gives a self- 
energy of order A^~". There are some corrections for finite to Eq. (15), 
but the most important higher-order ones (higher order in powers of 
compared to the mean-field) are corrections to the self-energy. The leading 
order correction to the mean-field involves diagrams in which a single link 
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Fig. 6. Average system resistance of SW networks with distance-dependent con- 
ductances as described in the text, for -p = 1.0 for different a values. The solid 
line segments are the measured slopes for a^^O, with their values shown in the leg- 
ends; these values should be compared to a/2, the exponent of the leading-order 
perturbative results [Eq. (19)]. 



appears twice; these diagrams involve summing over the length of the link 
and are mulitplied by the strength of the link squared: \i — For a;<l/2, 

this sum diverges for large and gives rise to corrections to the mean-field 
S, so that Eq = p(iV/2)~"/(l — a) is the mean-field value and the leading 
correction is E = p{N/2y/{l -a)- p[(Ar/2)-2"/(l - 2a)]/v^+ . . .. Using 
this correction to E, the corrections to the resistance are of order A^", and 
thus for a;<l/2 may be significant compared to the value in Eq. (19). The 
coefficient of these corrections becomes singular at q;=1/2 and for a>l/2, the 
self-energy E becomes non-local and acquires a momentum dependence, which 
may be shown to change the form of the higher-order corrections. 

Also note that this behavior is very different from that of the case where the 
strength (conductance) of the random links is uniform, but the probability of 
connecting two nodes, separated by a distance /, decays as [37,38,58]. There 
exists a finite region, 0<a<2, where the propagator, hence in the context of 
this paper, the average system resistance, remains finite in the limit of N^oo. 
In the present case, where the link-length distribution is uniform, but the link 
strength decays as the average system resistance is finite only for a=0. 

This contrasting behavior between the two different implementations of 

the random links (strength vs probability) is in accord with recent studies on 
phase transition on SW networks. Interacting systems often exhibit mean-field- 
like phase transitions [57,59,60,61,62,63,64,65], even for an arbitrarily small 
but nonzero density of random links added to a one-dimensional regular graph. 
However, in the case of the strength of the random links decaying in the above 
fashion, for the Ising model on SW networks, it was shown [66] that no 
phase transition occurs at any finite temperature for any nonzero a. 
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Summary and Outlook. — Wc obtained the scaling behavior of the effective 
resistance of SW networks. For uniform hnk conductances, we found that for 
an arbitrary small density of random links, the average system resistance is 
finite, and the two-point resistance, as function of the distance between the 
nodes, saturates exponentially fast to the same finite value. When the link 
conductance decays with the distance as /~", the average network resistance 
diverges with the number of nodes as N"/^. 

Ultimately, one is interested not only in the global transport or fiow charac- 
teristics of the network, but also in their effect on the local "components" , 
capacity limitations, and possible global network failures. In the context of 
resistor networks, the question of voltage landscapes in the network, or more 
specifically, the voltage-drop distribution across the links, can be addressed. 
Such a study can reveal the most vulnerable links/connections to be "blown" 
when increasing the overall load in the network. In particular, studying the 
properties of the extreme (largest) voltage-drops across the links in the net- 
work carries information on the weakest links of the network, and in turn, 
provides solutions from a system-design viewpoint. Fuse networks have been 
intensively studied on random percolating lattices with various applications 
to breakdown processes in condensed matter and materials systems, ranging 
from brittle fracture to dielectric breakdown [5,6,7,8,9]. Future work will ad- 
dress these questions from a general complex network vulnerability viewpoint 
[67,68,69]. 
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